Spin evolution of supermassive black holes and galactic nuclei 
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The spin angular momentum S of a supermassive black hole (SBH) precesses due to torques from 
orbiting stars, and the stellar orbits precess due to dragging of inertial frames by the spinning hole. 
We solve the coupled post-Newtonian equations describing the joint evolution of S and the stellar 
angular momenta Lj,j = 1...N in spherical, rotating nuclear star clusters. In the absence of 
gravitational interactions between the stars, two evolutionary modes are found: (1) nearly uniform 
precession of S about the total angular momentum vector of the system; (2) damped precession, 
leading, in less than one precessional period, to alignment of S with the angular momentum of the 
rotating cluster. Beyond a certain distance from the SBH, the time scale for angular momentum 
changes due to gravitational encounters between the stars is shorter than spin-orbit precession times. 
We present a model, based on the Ornstein-Uhlenbeck equation, for the stochastic evolution of star 
clusters due to gravitational encounters and use it to evaluate the evolution of S in nuclei where 
changes in the Lj are due to frame dragging close to the SBH and to encounters farther out. Long- 
term evolution in this case is well described as uniform precession of the SBH about the cluster's 
rotational axis, with an increasingly important stochastic contribution when SBH masses are small. 
Spin precessional periods are predicted to be strongly dependent on nuclear properties, but typical 
values are ~ 10 7 - 10 8 yr for low-mass SBHs in dense nuclei, ~ 10 8 - 10 10 yr for SBH masses 
~ 10 s Mq, and ~ 10 10 - 10 11 yr for the most massive SBHs. We compare the evolution of SBH 
spins in stellar nuclei to the case of torquing by an inclined, gaseous accretion disk. 

PACS numbers: Valid PACS appear here 
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INTRODUCTION 



An accretion disk fed by gas whose angular momen- 
tum is misaligned with that of the central supermassive 
black hole (SBH) will experience Lense-Thirring [1] pre- 
cession. Viscous torques near the SBH align the gas with 
the SBH equatorial plane [2]; farther out, the gas re- 
mains inclined, producing a constant torque that causes 
the SBH spin axis to precess. Such precession has been 
invoked as an explanation for changes in the direction of 
radio jets in active galaxies [3, 4]. Continued accretion 
of gas from a misaligned plane will eventually reorient 
the SBH, although the time required for realignment is 
uncertain [5]. 

Accretion disks are believed to be associated with only 
a small fraction of SBHs. Here we consider the more 
generic, and perhaps simpler, case of a rotating SBH 
embedded in a nuclear cluster of stars or stellar rem- 
nants. If the cluster has a net angular momentum that 
is misaligned with the SBH spin, a mutual torque will 
be exerted between stars and SBH, even if the spatial 
distribution of the stars is precisely spherical. 

In the simplest such model, the stars move indepen- 
dently of each other. Differential precession ( "phase mix- 
ing") will nevertheless cause stellar orbits near the SBH 
to distribute their angular momentum vectors Lj uni- 
formly about the spin S, decreasing the torque that they 



exert on the hole. The angular momentum associated 
with stars farther out can remain misaligned, leading to 
a forced precession of the SBH, similar to what occurs in 
the case of misaligned accretion disks. 

By solving the coupled post-Newtonian equations de- 
scribing a spinning SBH and a rotating cluster, we verify 
that such an outcome is possible, at least starting from 
certain initial conditions. However we find a second evo- 
lutionary mode as well, in which differential precession 
causes the inner system to reach alignment with the to- 
tal (spin plus orbital) angular momentum, resulting in a 
steady state with no subsequent precession of the hole. 

Stars also interact with each other gravitationally; 
these encounters lead to changes in stellar angular mo- 
menta, on time scales that can be short compared with 
Lense-Thirring times. Unlike changes due to frame- 
dragging, evolution of the Lj due to encounters is es- 
sentially random. There is a region near the SBH, the 
"sphere of rotational influence," in which encounter times 
are long compared with frame-dragging times. Within 
this region, stellar orbits precess uniformly, while out- 
side of it, changes in the Lj are due primarly to encoun- 
ters and are random. The size of this sphere varies from 
~ 1CP 3 pc in the nuclei of galaxies like the Milky Way, 
to ~ 10 1 pc in nuclei containing the most massive SBHs. 
We develop a stochastic model for the evolution of S 
that includes the effects of encounters on the Lj. In 
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this model, net alignment of the stellar angular momena 
with the SBH spin is less efficient, and the SBH typically 
continues to precess about the mean L of the stellar clus- 
ter, although its instantaneous precession rate can vary 
stochastically due to the stochastically changing Lj. 

Evolution of SBH spins to due torquing from stars has 
many parallels with evolution due to torquing from an 
accretion disk, surprisingly so given that one process is 
energy-conserving and the other is dissipative. We com- 
pare and contrast the two sorts of evolution in the "Dis- 
cussion" section, where we also summarize observational 
and theoretical evidence for nuclear rotation, and dis- 
cuss the implications of our results for the experimental 
determination of black hole spins. 

Throughout this paper we ignore the contribution of 
stellar captures to the evolution of S. 



II. SPIN-ORBIT EQUATIONS 



A Kerr black hole of mass M. has gravitational radius 



r g given by 



GM. H ( M, \ 

— - 48x10 (iwj pc (1) 



and spin angular momentum S, which we write in terms 
of the dimensionless spin parameter \ as 



S = X 



GM* 



lxl<i- 



(2) 



To lowest post-Newtonian (PN) order, the spin evolves 
(precesses) in response to torques from orbiting stars ac- 
cording to [6] 



dS 
~dt 



2G A w 



3=1 3 



(3) 



where rnj <C M, is the mass of the j'th star, whose in- 
stantaneous position and velocity relative to the SBH 
(assumed fixed at the origin) are {xj,Vj}, and tj = \xj\. 
Eq. (3) is invariant to the choice of spin supplementarity 
condition [6]. It can be written in the equivalent form 



S = u s xS, 



2G L 



where 



rrijixj x Vj) 



(4) 



(5) 



is the Newtonian angular momentum of the j'th star. 

We are mainly interested in changes that take place on 
time scales long compared with stellar orbital periods, P, 
where 



27ra 3 / 2 



VGMl 



2.96 



/ a V /2 f M. 



-1/2 



yr. (6) 



Here a is the orbital semimajor axis and mpc = 10~ 3 pc. 
Accordingly, each of the j terms on the right hand side of 
Eq. (4) can be averaged over the unperturbed (Keplerian) 
orbit, whose semimajor axis and eccentricity are a 3 and 



Using 



~=a- 3 (l-e 2 ) 



-3/2 



(7) 

and fixing Lj during the averaging, the spin evolution 
equation becomes 



S 



ljs x S 
2G 

c 2 

3 -j 



E 
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(8a) 
(8b) 



Henceforth averaging over the Keplerian motion will be 
understood unless otherwise indicated. 

Stellar orbits also precess in response to frame- 
dragging torques from the spinning SBH. Working again 
to lowest PN order and averaging over the unperturbed 
motion yields the standard expression for the Lense- 
Thirring [1] precession: 



Lj 



2GS 



c 2 a 3 (l-e 2 ) 3 / 2 



(9a) 
(9b) 



For fixed S, precession described by Eq. (9) has the 
form of uniform advance of the line of nodes, the latter 
defined as the intersection of the orbital plane with the 
equatorial plane of the SBH. We denote the nodal angle 
by fl; thus in the orbit-averaged approximation, flj = 
ujj. Orbits also experience precession of the argument 
of periastron due to both the Schwarzschild and Kerr 
components of the SBH metric, but such precession leaves 
the Lj unchanged. 

In the absence of interactions between stars, the cou- 
pled equations (8), (9) determine the joint evolution of 
the SBH spin and the stellar angular momenta. Con- 
served quantities include the total angular momentum of 
the system: 



J = S + '£L j = S- 



-'tot 



as well as 



S 

L 3> j 



1,...,A. 



(10) 



(11a) 
(lib) 



Neither S, L tot nor |L to t| is conserved. However conser- 
vation of S and J implies 



\J-S\<L tot < J + S. 



(12) 



Consider the case in which all stars have the same a 
and e; for instance, the orbits could lie in a circular ring. 
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There is no differential precession, and Eqs. (8), (9) can 
be written 



S = o; x S, 



uj x L t 



where 

U>LT 



2G 2 M 2 



c 3 a 3 (l - e 2 ) 3 / 2 
(7.0 x 10 5 yr) _1 



X 



(13) 

(14a) 
(14b) 



/ M. 
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In this special case, L to t is conserved, and both S and 
Ltot precess with frequency loo about the fixed vector 
J. The controlling parameter is = L to t/S. If <C 1, 
S ~ and L tot precesses about the fixed SBH spin vector 
at the Lense-Thirring rate; while if ^> 1, L to t ~ and 
S precesses about the fixed angular momentum vector of 
the stars with frequency 0wlt 3> wlt- 

This simple model might apply to the "clockwise stellar 
disk" at the center of the Milky Way, which has a mass 
~ IO^Mq, radius 0.04 pc < r < 0.5 pc, and mean orbital 
eccentricity - 0.2 [7-9]. Setting M,=4x lO 6 ^© [10], 
the implied is 



1/2 

~ o -1 ^CWD / -RcWd\ 



(15) 



consistent within the uncertainties with unity even if x 
is as large as 1. Evidently, the stars in this disk torque 
the SBH about as much as they are torqued by it. The 
mutual precession time is 



7T 



much longer than the ~ 10 7 yr age of the disk inferred 
from the properties of its stars, and also long compared 
with other physical processes that are likely to alter the 
stellar orbits (as discussed in more detail below). Never- 
theless, this example demonstrates that identified struc- 
tures near the Galactic center SBH can easily contain a 
net orbital angular momentum that exceeds S. 

The distribution of stars at distances < 0.1 pc from 
Sgr A* is poorly constrained [11-13], but the total stel- 
lar mass in this region is almost certainly large compared 
with the ~ 10 4 7W© associated with the clockwise disk. 
Given the strong (oc r~ 3 ) radial dependence of the frame- 
dragging torques, even a modest degree of net circulation 
of the stars in this region could therefore induce a pre- 
cession of the SBH on time scales very short compared 
with the time of Eq. (16). 

We emphasize that there is no need for the torquing 
stars to lie in a geometrically flattened structure: accord- 
ing to Eqs. (8), all that is needed is a non-random orien- 
tation of the orbital angular momentum vectors, which 



occurs even in a precisely spherical nucleus if there is a 
preferred sense of orbital circulation. 

In general, different stars will have different a and 
e, implying different rates of nodal precession. Close 
enough to the SBH, orbital precession times will be short 
compared with the precessional period of the SBH, and 
the orbits will tend to distribute their angular momen- 
tum vectors uniformly about the instantaneous S. The 
net torque from these stars will then fall essentially to 
zero, and continued precession of the SBH will be driven 
by stars farther out. We expect the radius separating 
stars in these two regions to be roughly the radius con- 
taining a total stellar angular momentum equal to S. We 
estimate that radius in the following section, after first 
presenting observationally motivated models for stellar 
nuclei. 



III. SPHERICAL NUCLEI 

Most of the distributed mass at distances r < 0.1 pc 
from the Milky Way SBH is believed to be in the form 
of stars much older than the stars in the clockwise disk. 
The spatial distribution of these stars is believed to be 
approximately spherical [14], with at least a modest de- 
gree of circulation [12, 15]. 

A simple model for the distribution of mass near the 
center of a spherical galaxy is 



p(r) 



r 

Po \ — 

Jo 



(17) 



Near the SBH (but not so near that relativistic correc- 
tions are required) , the gravitational potential is 



$(r) = - 



GM. 



(18) 



and orbits can be characterized by their semimajor axes 
and eccentricities, as in the previous section. If the stellar 
velocity distribution is assumed to be isotropic and sta- 
tionary, and if stars are distributed along orbits uniformly 
with respect to mean anomaly, the joint distribution of 
a and e that generates the density (17) is 



N(a,e) dade — N a 2 7 daede. 



(19) 



The relation between No, po and ro is easily shown to be 

7 > 1/2 (20) 



8^ 3 / 2 r( 7 + i) 7 



2t r(7-l/2)' 

where m* is the stellar mass, assumed the same for all 
stars. Values of 7 less than 1/2 are not achievable if the 
velocity distribution is isotropic [16]; we do not consider 
that possibility here, and in the modelling that follows, 
7 will be restricted to the range 1/2 < 7 < 3. 

The relations (18) - (20) are valid at radii smaller than 
the SBH influence radius r m , customarily defined as the 
radius enclosing a stellar mass equal to 2M # : 



M*(r < r m ) = 2M.. 



(21) 
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A spherical cluster will exhibit net rotation if un- 
equal numbers of stars (at each a and e, say) circulate 
in a clockwise vs. counter-clockwise sense about some 
axis. For instance, if one-half of the orbits in a spher- 
ical cluster with initially isotropically-distributed veloc- 
ities have their velocity vectors reversed such that all 
angular momentum vectors point toward the same half- 
sphere, the total angular momentum of the ensemble will 
be |.Ltot| = | S I A? I- Henceforth we characterize the net 
rotation of a spherical cluster by the factor /, defined 
as the fraction of orbits that have been "flipped" in this 
way; < / < 1/2 and / is assumed to be independent of 
a and e. 

Characterizing the rotation in this way is "conserva- 
tive," in the sense that a geometrically flattened nuclear 
cluster (e.g. a disk), or a spherical cluster consisting of 
only circular orbits (an "Einstein cluster" [17]), can have 
a larger net angular momentum for the same radial dis- 
tribution of mass. 
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FIG. 1: Schematic representations of the two nuclear mod- 
els considered here, (a) Low-mass galaxy with nuclear star 
cluster (NSC), (b) High-mass galaxy with core. 

Observed galaxies appear to fall into one of two classes 
in terms of the parameters that define their stellar dis- 
tribution at r < r m [18]. Massive spheroids - elliptical 
galaxies, or the bulges of spiral galaxies - with total lu- 
minosities greater than ~ 10 10 - 5 L© have "cores," regions 
of size ~ r m where the stellar density rises slowly toward 
the SBH. In these galaxies, the observed, mean relation 
between M. and r m is approximately [19] 

rm ^ 35 (los^J pc ' a " - 56 (22) 

and the index 7 that defines the central density increase 
varies from <~ at the highest luminosities to <~ 2 or 
<~ 2.5 at the low- luminosity end of the range, albeit with 
substantial scatter [20, 21]. The SBHs in these galaxies 
have masses 10 7 - 5 M @ < M, < 1O 9 5 7W . 

Less luminous spheroids often exhibit dense cen- 
tral mass concentrations, called "nuclear star clusters" 
(NSCs). The sizes of NSCs are also comparable with r m 
(assuming that the host galaxies contain SBHs), although 
these structures are too compact to be well resolved in 
galaxies beyond the Local Group. The best-studied case 



is the Milky Way, in which the stellar density appears to 
follow p(r) ~ r" 1 ' 8 inside ~ 5 pc, compared with a SBH 
influence radius of ~ 2.5 pc [14, 22]. The high densities of 
NSCs imply short time scales for equipartition of orbital 
energies [23], and one expects the densest NSCs to ex- 
hibit mass segregation, i.e. the heavier bodies should be 
more strongly concentrated toward the center than the 
lighter bodies. The heaviest bodies are expected to be 
stellar-mass black holes (BHs), the end products of stars 
with initial masses m* > 30 Mq whose main sequence 
evolution requires only a few million years; BH masses 
are believed to be in the range 5.M© < < 2QAi Q 
[24], compared with a main-sequence turnoff mass of 
~ 1M.Q. When energy equipartition is satisfied, the 
lighter population is predicted to follow p(r) ~ r - 3 / 2 
at r < 0.2r m while the BHs obey the steeper relation 
p <~ r~ 2 [25, 26]. Detailed dynamical models of the 
Galactic center [27, 28] suggest that if the nucleus is older 
than an energy equipartition time, about one-half of the 
distributed mass inside 0.01 pc would be in the form of 
main-sequence stars and one-half in BHs, with a smaller 
mass fraction in neutron stars and white dwarves. How- 
ever it is currently unclear whether the Milky Way NSC 
has a relaxation time short enough for gravitational en- 
counters to have produced such a distribution in 10 Gyr 
[13] and the distribution of observed giant stars (with 
masses ~ 1 — 3.M0) is much flatter than predicted in the 
relaxed models inside ~ 0.5 pc [9, 29, 30]. 

In what follows, the central regions of bright and faint 
galaxies will be parametrized in different ways. Nuclei 
of bright galaxies, with M, > 10 7 - 5 A4q, are assumed to 
follow Eq. (17) at r < r m , with r m determined by M. 
via Eq. (22). The distributed mass interior to r in these 
galaxies can be written 

(\ 3—7 
— J (23) 

/3 = 1 -a(3- 7) w -O.68 + O.567. 

Mass segregation is expected to be unimportant in the 
nuclei of giant galaxies so we set m* = a typical 

value for an old stellar population. 

In the case of galaxies with M, < lO 7 ' 5 ^©, the distri- 
bution of mass at r < r m is less certain. We parametrize 
these nuclei in terms of both M. and M .i, the latter 
defined as the mass in stars or stellar remnants inside 
r = 0.1 pc. If the power-law dependence of density on 
radius in these galaxies were to extend outward as far as 
r m , and if r m varied with M. as in bright galaxies, then 

Mcr =2M.(^) 7 " 3 (24a) 

~ 2x1 ° 3+, MtwU ■ (24b) 
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Eq. (24b) could be taken as a rough guide to the expected 
value of M i, but both M .i and 7 will be considered free 
parameters. We expect 1 < 7 < 2 for these nuclei; the 
stellar mass will be set cither to IMq (stars) or IOMq 
(stellar BHs). 

In both kinds of nuclei, rotation will be parametrized 
in terms of the fraction of flipped orbits, /, defined above. 



For low-luminosity galaxies, we find 




1 000 10' 10' 

M, (r<0.1pc) (M e ) 

FIG. 2: Characteristic radii as a function of nuclear param- 
eters. aL, Eqs. (27)-(28), is the semimajor axis containing 
a total angular momentum equal to S, computed assuming 
/ = 1/2 and \ = I (maximum rotation of SBH and stellar 
cluster), ax, Eqs. (55) - (56), is the radius of rotational in- 
fluence of the SBH, assuming \ = 1- Fhe left panel assumes 
M. = 10 6 Mq and r g = 4.8 x 1CT 8 pc; the stellar density is 
parametrized in terms of M0.1, the mass within 0.1 pc, and 
7, the power-law index; solid lines are for 7 = 1 and dashed 
lines for 7 = 2. In the case of <zk, two values are assumed for 
the stellar mass: m* = IAAq (thin lines) and m* = IOAAq 
(thick lines). Fhe right panel, for massive galaxies, assumes 
the relation (22) between M. and the influence radius r m ; 
solid, dashed and dotted lines are for 7 = 5/8, 1 and 3/2 
respectively. Fhe curves for hk in the right panel assume 
m* = IMq. Fhe radius of tidal disruption of a solar- type 
star falls below the lower boundary of both panels. 



Fhe total angular momentum associated with stars 
whose semimajor axes are less than a is 

L tot (a) = {fe L )Y,m j [GM.a j {l-e])] 1 ' 2 (25a) 

a,j <a 

-+ (fe L )(GM.) 1/2 N a m* f daa 5 ' 2 ^ 

Jo 

f 1 1 
x / dee(l-e 2 ) 

Jo 



1/2 



(25b) 



3(7 -2 7 ) 



(fe L )(GM.) 1/z N m*a 



V 2 JV„ m .„7/2- 



(25c) 



with eL a unit vector in the direction of L to t ■ We define 
aL such that 



Aot(aL) = i> = x - ■ 



(26) 



0.1 pc 



7/2-7 



1.5 x 10 



_ 2 x2^(7-27)r( 7 -l/2) 



/ M .! \ V M. 

X \10 6 M G 
while for bright galaxies, 
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3/2 
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(27) 



35 pc 



7/2-7 



3.9 x 10~ 5 



X 2^(7 - 2 7 ) r( 7 - 1/2) 
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(- 



M. 



\IQ S M( 



1/2+0(3-7) 



r(7 + i) 

(28) 



Figure 2 plots ol as a function of nuclear parameters. In 
massive galaxies, and for x/f ~ 1 ; 



10" 2 r n 



Fhe approximate radius of tidal disruption of a Solar- 
type star is [16] 



n 



that is 



9-8 X ID"' (-r^- 

V1O 8 M 



1/3 



mpc 



(29) 



V1O 8 M 



-2/3 



This radius is small compared with all radii relevant to 
the spin evolution of SBHs. Compact remnants would 
not be affected by tides from the SBH at any radius 
greater than r g . 

Based on the arguments in the preceding section, we 
expect stars at r < to precess about the SBH in a time 
short compared with the precession time of the SBH. 



IV. SPIN-ORBIT EVOLUTION 

The focus in this section is on the large-iV, or "col- 
lisionless," limit, appropriate for giant galaxies in which 
the central density is low and time scales for gravitational 
interactions between stars are long. (A more precise cri- 
terion is given in §V.) Accordingly, the number of stars in 
the numerical integrations was chosen to be large enough, 
typically N = 10 6 , that discreteness effects were small; 
otherwise the value of N is unimportant. 

Assuming a density law (17), the coupled evolution 
equations (8) and (9) admit of straightforward scaling 
relations. If the distributions of orbital eccentricities 
and inclinations are invariant under the rescaling, we can 
write 



LOj cx M 2 aJ 3 X . 



-(T +1 /2)do, (30a) 
(30b) 
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Consider first the case r m oc a w 0.56 that was 

adopted for luminous galaxies. Setting r = r m in Eq. 
(17) gives p = p(r m ) oc M.r m 3 ex M. 1 " 3 ". Then 



,,(3-5a)/2 

w s oc Mi 
Wj oc M. 2 - 3a X . 

Scaling M. and x independently as 

M. ->CiM., X^C* 2 x 

then yields 

£y(3-5a)/2 



(31a) 
(31b) 

(32) 

(33a) 
(33b) 

(34) 

if the unit of time, [T], is to scale the same way in both 
evolution equations. With this choice, 



u s oc 



UJj 



oc C?- 3q C 2 . 



Evidently we require 



C 2 = C? 



i(a-l) 



Cr 



0.2 



[T] oc Cf a - 3)/2 ~ C? 



(35) 



since a w 0.56 w 3/5. 

In the case of low-luminosity galaxies, the nuclear den- 
sity was specified by the independent parameter Mo.i , the 
stellar mass inside 0.1 pc. Defining a third scale factor 
as 



o.i 



C 3 M . 



it is clear that 

,-,1/2 ,-,-1 

Wj oc C1C2 
and a common unit of time requires 



C 2 



c\ 3/2 C3- 1 . 



(36) 

(37a) 
(37b) 

(38) 



For both sorts of rescaling, the condition % < 1 implies 
limits on the values of C\ and C3. 

Integrations of the coupled equations (8), (9) were car- 
ried out using a 4(5) order Runge-Kungc routine with 
adaptive time steps [31]. Monte-Carlo initial conditions 
for N stars were first generated from Eq. (19) assuming 
a random distribution of orbital planes, i.e. an isotropic 
velocity distribution. An upper limit, a max , was imposed 
on a, and a lower limit, r Pim i n , on the radius of orbital 
periapsis r p = a(l — e). A fraction / of the orbits at each 
(a, e) were then "flipped" (the sign of Lj was changed) in 
order to give the cluster a net rotation about the z-axis. 

For these initial models, the spin precession vector u>s, 
Eq. (8b), is given by 



3 

2G 3 / 2 M. 1/2 



1 1/2 



of (1 - e 2 ) 3 / 2 



(39a) 



-(/e L )iY m* 



da ede 



aV2+7(l 




time (yr) 

FIG. 3: Evolution of the SBH spin precession rate in a set of 
integrations with r Pim i n = (0.03,0.1,0.3) mpc and a max = 30 
mpc. The other model parameters are specified in the text. 



where m 4 is the mass of one star and is a unit vector 
in the direction of L tot . The integral (39b) diverges as 
the integration limit in a tends to zero for 7 > 1/2, or 
as the limit in e tends to one. A lower limit could be 
placed on a(l — e) by the requirement that stars come 
only so close to the SBH before being captured or tidally 
disrupted. But as noted above, one expects the net an- 
gular momentum of stars at small radii to align quickly 
(on a time scale much shorter than the time for changes 
in S) with S, reducing their contribution to dS/dt. 

That this does indeed occur is illustrated in Figure 3, 
which shows integrations of a set of models that differ 
only in the choice of r p . m i n . The models have M. = 
1O 6 A4 , x = 1, a max = 30 mpc, / = 1/2, 7 = 1, and N = 
10 6 . The SBH spin axis was oriented initially at an angle 
of 60° with respect to L tot . For these parameters, <zl ~ 5 
mpc and wlt(ol) ~ 1 x 10~ 8 yr -1 . The initial conditions 
with smaller r p , m i n have larger initial ljs- However the 
torque from the inner stars decays on a time scale of 
order the Lense-Thirring time for the innermost orbits 
as their angular momentum vectors distribute themselves 
uniformly about S, and S hardly changes in this time. 

The long-term evolution of the models in Figure 3 con- 
sists of precession of the SBH about J w L to t ■ It turns 
out that a second evolutionary mode is possible in spheri- 
cal models like these. This is illustrated in Figures 4 and 



5, based on a cluster with parameters M, = 10 M.q 



X = 1, 7 = 1, / = 1/2, a n 



100 mpc, r % 



1 



(39b) 



mpc, and a total stellar mass of 1O 5 A^0. For this model, 
Mo.i w 6 x 1O 4 A40 and w 15 mpc. The integrations 
shown in Figures 4 and 5 are from a sequence in which 
the initial angle, 8 , between S and L to t was varied in 
steps of 10°, from 10° to 170°. For 8 > 45°, evolution 
at late times consists of nearly uniform precession of the 
SBH spin axis about J, as in the integrations of Figure 3. 
However if O < 45°, precession continues only for a sin- 
gle cycle or less, after which the vectors S, L tot and u>s 
are nearly aligned and precession essentially stops. 
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FIG. 4: Evolution of S, -Ltot, and u>s in a cluster where the initial SBH spin axis was offset by Oo = 40° from the stellar 
angular momentum vector. The other parameters of the model are given in the text. In the upper panels, the open/filled 
circles indicate initial/final times respectively. This is an example of damped precession: the vectors S, Ltot and u>s reach a 
common orientation after roughly one precession cycle. Qualitatively the same sort of evolution occurs for < 8o < 45°. 
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other parameters of the initial models are given in the text. 



is, have very little scatter in their individual precession 
frequencies u)j , steady precession without alignment may 
persist even for L to t <C S. 

V. INFLUENCE OF GRAVITATIONAL 
ENCOUNTERS ON SBH SPIN 

Times associated with Lense-Thirring precession about 
a SBH arc long, and over such long time scales, stellar 
orbits can evolve in response to other influences. Here 
we consider how (Newtonian) gravitational interactions 
between stars would alter the evolution of SBH spins. 
These interactions are expected to be most important in 
the dense nuclei of low-luminosity spheroids; we derive 
more exact criteria below. We continue to assume that 
S evolves according to Eqs. (8), but we now allow for the 
possibility of other terms in the the evolution equations 
for the Lj. 



Evolution of the second sort, or "damped precession," 
which leads to almost complete alignment of SBH spin 
with L t ot , is not excluded by the conservation laws (10), 
(11), and in principle could occur for any initial condi- 
tions. In practice, we found that it occurs only when 9q 
is sufficiently small. The critical angle, # cr it, separating 
the two evolutionary modes was found to depend on the 
other model parameters. Figure 6 shows the dependence 
of # cr it on the mass of the stellar cluster, when the other 
initial parameters are the same as in Figures 4 and 5. 

A large number of such integrations revealed that the 
two modes of evolution illustrated in Figures 4 and 5 are 
generic. Roughly speaking, the system may end up in 
one of two distinct states: 

• Aligned S, L to t and J; 

• Uniform precession of both S and L to t about a 
fixed axis, essentially the axis of total angular mo- 
mentum J = S + L to t- 

In the latter case, typically the angle 9 between S and 
Ltot decreases from its initial value 6> , but settles at 
some non-zero average value after a couple of precessional 
periods. As noted above, the overall precession frequency 
may be estimated as the Lense-Thirring time for stars at 
the radius such that the total angular momentum of stars 
within this radius is equal to |S|. This frequency depends 
only weakly on the angle 9o provided that 9o 7^ 0. 

If the total angular momentum of the stars is less than 
S, essentially all the orbital Lj end up aligned or counter- 
aligned with the SBH spin, depending on whether 9o is 
greater or less than tt/2. This situation is unlikely to 
be relevant for galactic nuclei, since there will always be 
enough stars sufficiently far from SBH that their total 
angular momentum exceeds S, although the precessional 
times associated with distant stars may be long. In ad- 
dition, if the stellar orbits are initially concentrated in 
a small interval of radii on nearly-circular orbits, that 



A. Encounter time scales 

To a first approximation, the force from the stars 
can be modelled by approximating their distribution as 
spherically symmetric and stationary. The addition of a 
spherical component to the otherwise Keplerian poten- 
tial of the SBH results in an advance of orbital periapsis 
of each star ("apsidal precession") at an orbit-averaged 
rate given by [16] 



vm = -v r G M (e,j) (l - e 2 ) 



2\l/2 



M»(a) 
M. 



with associated time scale 



tM = 



7T 



p 



M. 

M*(a) 



(40) 



(41) 



Here, v r = 2ir/P is the Keplerian (radial) frequency, 
M+(a) is the mass in stars within radius r — a, and 
Gm ~ 1 is a weak function of 7 and e [16]. Adopting our 
parametrization for low-mass galaxies, this becomes 



1.5 x 10 5 / M. 



(l-e 2 ) 1/2 \Af .i 
while for high-mass galaxies, 



O.lpc 



7-3/2 



yr (42) 



2.3 x 10 4 



(l-e 2 ) 1/2 V!PC 



3/2 



7-3/2 



yr. (43) 



This "mass precession" leaves the orbital plane, and 
hence Lj, unchanged and so does not directly affect the 
evolution of S as given by Eq. (8). The same is true 
for the in-plane precession due to the Schwarzschild and 
Kerr parts of the SBH metric; in the orbit-averaged, post- 
Newtonian approximation, the time scale asociated with 
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the former precession, which always dominates the the 
Kerr contribution, is 



*s = 



7T (1 - e 2 )a 5 / 2 c 2 
3 (GM.) 3/2 



1.0 x 10 9 (l - e 2 ) 



5/2 



O.lpcy ^1O 6 M 



-3/2 



yr. 

(44) 



This "Schwarzschild precession" is more rapid than 
mass precession when 



(1 - e 2 ) 3/2 ( — ) <1 



where 



a s M*(r < a s ) = 3M.r g . 
For low-mass galaxies this is 



,0.1pc/ Mo.iO.lpc 
and for high-mass galaxies, 

4—7 



3 Tg 

2r„ 



(45) 



(46) 



(47) 



(48) 



For example, setting 7 = 2 in the first relation gives 



a s « 1.2 



/ M» \ ( Mq.i 



'1/2 



mpc. 



(49) 



While not directly affecting the Lj, these two sources of 
precession are important in setting the time scale for ran- 
dom fluctuations in the orbital eccentricities, as discussed 
in more detail below. 

Newtonian perturbations can also mimic frame- 
dragging by changing the orientation of orbital planes. If 
such changes occur on a time scale that is short compared 
with the Lense-Thirring precessional time, the evolution 
of orbital orientations will be determined essentially by 
the Newtonian perturbations [32]. We expect this to 
be the case for stars that are sufficiently far from the 
SBH, since frame-dragging time scales increase rapidly 
with distance (Eq. 14b). 

Here we focus on a generic source of non-spherically- 
symmetric perturbations: resonant relaxation (RR), the 
changes in L that result from the finite- N asymmetries in 
an otherwise spherical cluster around a SBH [33]. (Other 
possible sources of non-sphericity, ignored here, include 
a large-scale distortion of the nuclear potential or "bar" 
[34], or a distant massive perturber [35].) In what follows, 
we call the evolution of orbital planes due to these mutual 
torques "2d resonant relaxation," or 2dRR [94]. 



Under 2dRR, orbital orientations change in a charac- 
teristic time [16] 



?2dRR 



P M. 1 

4.7 x 10' 
M. 



(50) 



10 6 



a 
mpc 

' N 

To 2 



3/2 



/ M. 

\lO e M Q 



-1/2 



-1/2 



yr 



where P = P{a) = 2-K/v r is the radial (Kepler) period 
and N = N(a) is the number of stars at r < a. In a time 
~ Tmrr, orbital planes will have essentially randomized 
due to the mutual torques [95]. 

The condition that frame dragging causes orbital 
planes to precess more rapidly than they are changed 
by the mutual torques is 



(x>LT 



or equivalently [96] 
(1-e 2 ) 3 



< 



16x 2 / M, 
N(a) V^T 



(51) 



(52) 



Orbits satisfying this condition will be said to be in the 
"collisionless" regime: to a first approximation, their an- 
gular momenta evolve in accordance with Eq. (9), unaf- 
fected by perturbations from other stars. 

The condition (52) can be expressed in terms of a char- 
acteristic semimajor axis, cik, as 

)6— 7 
< 1. (53) 

We call aK the "rotational influence radius" of the SBH. 
To solve for ok, we write N(a) for each of the two types 
of nuclear model defined in 5111 as 



(1-e 2 ) 3 



N(a) 



N(a) 



Mi 



0.1 



0.1 pc 

3-7 



3-7 



M. < 10 7 b M, (54a) 



M, > 10 A^ Q . (54b) 



These approximate expressions are adequate given the 
approximate nature of Eq. (52). In the case of low-mass 
galaxies, Eqs. (52) - (54a) yield 



ax 
0.1 pc 

/ M. 



6 — 7 



1.8 x 10~ n x 2 



(55) 



VlO 6 7W y \W 5 M Q J \IM Q/ 

while for high- mass galaxies, Eqs. (22), (52) - (53) and 
(54b) give 

6—7 

(56) 



35 pc 



2.1 x 10~ 12 x 2 



/ M. 



4+a(3-7) 



10 



with a w 0.56. Figure 2 plots ok as a function of nu- 
clear parameters. In low- mass galaxies, «k C «l; as M, 
increases, ok can approach cil- In the latter case, we 
expect the net angular momentum associated with stars 
inside the rotational influence sphere to be comparable 
with 5. 
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FIG. 7: 9k is the angular momentum associated with stars 
in the collisionless regime, expressed as a fraction of the SBH 
spin S = xGMljc. Fhe left panel assumes M. = 10 6 Mq, 
X = 1 and / = 1/2 (maximal rotation of SBH and stellar 
cluster); the stellar density is parametrized in terms of A/o.i, 
the mass within 0.1 pc, and 7, the power-law index; solid 
lines are for 7=1 and dashed lines for 7 = 2. Fwo values 
are assumed for the stellar mass: m* = IMq (thin lines) and 
m* = IOAAq (thick lines). The right panel, for massive galax- 
ies, assumes the relation (22) between M. and the influence 
radius r m ; solid, dashed and dotted lines are for 7 = 5/8, 
1 and 3/2 respectively. Thick lines assume X = 1, / = 0.5 
(rapidly rotating cluster) and thin lines assume x = 1, / = 0.2 
(slowly rotating cluster) . All curves in the right panel assume 
m* = IMq. 



Define 



Lk 

s 



(57) 



where Lk is the angular momentum associated with stars 
that satisfy (53). We compute Lk from Eq. (25b) after 
modifying the integration limits to respect the condition 
(53). The result is 



e. 



for bright galaxies, and 



for faint galaxies, where 
H{>y) = 

8^ (3-7) r( 7 + i) 

3 27 (7 - 1) r( 7 - 1/2) 
for 7 7^ 1, and 




3-7 




3-7 



(58) 



(59) 



6 — 7 
7-2 7 



a K 



(60) 

(l-7)/2 



F «=15 + 3 l0g 



OK 



(An upper cutoff to a is only required when 7 < 1 due 
to a weak divergence of the integral; a natural choice is 
flmax = r m since the expressions for N(a, e), etc. are only 
valid at r < r m .) 

Figure 7 plots 0k as a function of nuclear parameters. 
As expected, for high-mass galaxies, and for / ss 1/2, Ok 
is of order unity, scaling as fx ^ for fixed M.. As 

M, is decreased, 6k falls as well, although it can still be 
appreciable, 0.01 < 6k ^0.1, in low-mass galaxies with 
dense nuclei, 7 w 2. 

Figure 7 suggests that in galaxies with large M., the 
joint evolution of S and Lj will be similar to the evolu- 
tion described in the previous section, in the sense that 
mutual stellar interactions can be neglected. As M. is 
decreased, the angular momentum associated with stars 
in the collisionless regime drops compared with S. In nu- 
clei with sufficiently small M, , most of the torque acting 
on the SBH is likely to originate in stars whose orbits re- 
spond to each other on a shorter time scale than the local 
Lcnsc-Thirring time. As a result, the angular momentum 
vectors of these stars will be unable to align around S 
as in the collisionless case. We argue in the next sec- 
tion that the result can be substantially higher rates of 
sustained SBH precession. 



B. Stochastic model for the evolution of u>s 

In principle, the combined effects of gravitational self- 
interactions and spin-orbit torques could be directly sim- 
ulated using an A-body algorithm [32]. However the ra- 
tio between Kerr precessional times and orbital periods is 
so great that such direct simulation would be expensive 
for any reasonable N. 

An alternative approach would be to incorporate the 
effects of star-star interactions by modeling the evolu- 
tion of each of the Lj as a random walk [36-38]. How- 
ever, interactions between stars must conserve L to t, as 
well as being constrained in less obvious ways by the fact 
that the torques are mutual. Approximating the evo- 
lution of each star's angular momentum as an indepen- 
dent stochastic process, independent of the changes in 
the other Lj, would fail to capture these essential con- 
straints. 

Since the effects of the Lj on S appear only through 
u>s, and since the time scales for changes in the Lj due to 
self-interactions are typically short compared with spin- 
orbit time scales, it is reasonable to separate the prob- 
lem into two parts: asking first how u>$ varies as the 
stars interact with one another, ignoring the effects of 
spin-orbit torques; then using this knowledge to predict 
how S would evolve in response to the fluctuating us- 
We first explore this model, then present a more careful 
justification below. 

Figure 8, based on a direct integration of the A^-body 
equations of motion for 100 point masses (stars) orbit- 
ing about a massive particle (SBH), illustrates how u>s 
evolves due to star-star interactions. The integrator 
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FIG. 8: Evolution of u>s in the A^-body integration described 
in the text. Units of u> are inverse years. 
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FIG. 9: Distribution of the u>s values measured in the N- 
body simulation of Figure 8. Solid lines are Gaussian fits 
with a x , v , z = {5.8,7.0,4.8} x KT 10 vr" 1 . 



[39, 40] included 1PN terms in SBH-star interactions; 
spin-orbit terms were omitted. Initial conditions were 
generated according to Eq. (19), with M, — lO 6 ^©, 
to* = IOAI0, a max = 10 mpc, r Pim i n = 1 mpc, 7 = 1 and 
/ = 1/2. Rotation of the cluster was initially about the 
z-axis. 

Total angular momentum, L tot (or rather, its 1PN ana- 
log [41]) is conserved in these A^-body integrations. The 
spin precessional vector is not conserved; but since ujs is 
a weighted sum of the Lj , and since ^2 Lj is conserved, 
exchange of angular momentum between stars tends on 
average to leave unchanged. However, Figure 8 shows 
that each component of fluctuates about its mean 
value in an apparently random fashion. The amplitude 
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FIG. 10: Circles are auto-correlation functions of the time 
series plotted in Figure 8; horizontal axis is the lag time. 
Arrows mark the computed correlation times. Solid lines are 
Eq. (62). 



of these fluctuations is approximately constant over time, 
giving each time series the appearance of a stationary 
stochastic process [42]. 

Assuming stationarity, it is reasonable to calculate the 
distribution function of the fluctuations at any given time 
by binning together the events from all times. The results 
are shown in Figure 9 where the distributions have been 
fit to Gaussian functions. 

The time scale associated with stochastic fluctuations 
in in the A^-body integrations can be found by com- 
puting the autocorrelation functions (ACF), defined as 



Ri 



ws,,] [us,i(t' +t) - ^sA dt' 



Sol u s. 



(61) 

Here, uis,i is the ith. component of u)s, ^>s,i is its time- 
averaged value, and < t < T is the elapsed time in the 
A^-body integration. Figure 10 shows that the measured 
ACF's are reasonably well fit by exponential functions: 



R u {t) wexp(-t/Ti), i — {1,2, 3}, 



(62) 



with Tj « 1.5 x 10 5 yr. 

One expects the autocorrelation time for u)s{t) to be 
similar to the characteristic time associated with changes 
in the Lj. One such time is the 2d resonant relaxation 
time defined in Eq. (50). There will also be variations 
in ujg due to changes in orbital eccentricities. In the so- 
called "coherent RR" regime, defined as At < {tM,ts}, 
changes in e occur in a characteristic time ~ T-mrr; while 
in the "incoherent RR" regime, i.e. At > {£m,£s}, the 
associated time is longer than T 2( jrr (Appendix). Hence 
?2dRR is the shortest of the variability time scales and we 
can safely associate r with it. For the A^-body models, 
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Eq. (50) states 



again as 



T 2 dRR « 5 x 10 4 



1 mpc 



(63) 



quite consistent with the autocorrelation times measured 
in the TV-body simulation given that lmpc < a < 
10 mpc. The longer time scale associated with random- 
ization of orbital eccentricities is given by Eq. (A3): 



T R R jM (a) « 3 x 10 6 



1 mpc 



\ 3/2 



(64) 



for a > as ~ 0.00 mpc. 

The variance in the components of ws, 



° % = [us,i(t) - us,i] 2 



(65) 



can also be estimated given the known properties of the 
initial model. Begin by rewriting us, Eq. (8b), as 



us,i{t) 



3 

2G 3 / 2 M. 1/2 m 3 
Aj,i(t) = [uLj(t)-«i] 



(66a) 

(66b) 
(66c) 



with UL,j{t) a unit vector in the direction of Lj and Ui 
a unit vector in the direction of the zth coordinate axis. 
In general, each of the variables {e-j.Xj^} will change 
stochastically due to star-star interactions, and all of 
these changes will contribute to the variance of u>s- A 
lower limit on that variance follows from assuming that 
resonant relaxation causes changes only in the orbital 
planes and that the ej are approximately constant. Then 



a l = var [ui s ,i 



E ( 1 x var t A ^i 

3 
3 

AG 3 M.m 2 



E 



(67a) 
(67b) 

(67c) 



assuming var [A-,] w 1. For a cluster containing orbits 
with a single (a, e), the right hand side is ~ w| max /A, 
where ws. max is the magnitude of in a maximally- 
rotating cluster with / = 1/2. Then a w ws.max/v^/V- 

According to Eq. (64) , orbital eccentricities should also 
change substantially over the integration period, partic- 
ularly for orbits of small a. The variance in wg should 
therefore contain a substantial contribution from changes 
in the ej , and in fact the formula just derived underpre- 
dicts the variances observed in the A-body simulations 
(Figure 9) by a factor of a few. We estimate a allowing 
for changes in the ej as follows. Rewrite Eq. (8b) yet 



A3 

X 3\i 



2GV 2 M l J 2 mi 



2 5/2 
Ca 3 



l-e 2 



(68a) 

(68b) 
(68c) 



where both \j and ej are allowed to be functions of time. 
Assuming uncorrelated changes, the variance in us,i is 



3 

3 



A 2 var [Xj] , 



(69a) 



var [A,] = (A j;i ) var [(1 - e 2 )' 1 } 
+ (l- ej 2 )-ivar[A^] 
+ var [(1 - e 2 )- 1 ] var [A j;i ] . 



(69b) 



We estimate the quantities on the right hand side of 
Eq. (69b) by assuming that resonant relaxation main- 
tains a "thermal" distribution of eccentricities at every a 
[37], i.e. that 



N(e)de w 2ede 
for < e < e max (a), e max < 1. Then 



(1 



var 



[ln(l-eL x )] 2 - 



(70) 



(71) 



(72) 



We identify e max (a) with 1— r Pjm ; n /a. We likewise assume 
that orbital planes are randomized, so that varfAj^] = 
1/3, and Xj i = f in the case that ttj || L tot and zero 
otherwise. Finally, since e max w 1, it is reasonable to 
ignore the logarithmic terms, yielding 



2G 3 M. 



p.mm 



p,min 



(73) 



for i = 1,2,3. 

Applied to the A-body models in Figure 9, Eq. (73) 
yields Oi k> 6 x 10~ 10 yr -1 , in good agreement with the 
values obtained via the Gaussian fits to the A-body data. 

As long as the characteristic time for changes in ec- 
centricity is shorter than the other times of interest, Eq. 
(73) is the appropriate expression to use for a. This will 
turn out always to be the case in the examples presented 
below. 

A theorem [42] states that a stationary random pro- 
cess with a Gaussian probability function and an expo- 
nentially decaying autocorrelation function is necessarily 
an Ornstein-Uhlenbeck [43] process. The latter is defined 
as having a transition probability between two states, y\ 
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FIG. 11: Ornstein-Uhlenbeck realization of ws(t) using (r, cr) 
derived from the iV-body simulations in Figure 8. 



and j/2 (given here by two values of 0)5)7 at times t\ and 
ti that obeys 



T{yi\y\) 



V27r(l- e -2A) 



cxp 



(2/2 - J/ie- A ) 2 



2(l-e- 2A ) 



(74) 

where A = {ti — ti)/r and r is defined as in Eq. (62). An 
OU process X(t) with mean value X can also be defined 
as the solution of the Langevin equation, 



X = - 7 [X(t)-X] +AT(t) 



(75) 



if 7 = r 1 and if Af(t) is a Gaussian random variable 
having the properties 



M(t) = 0, 
N{t)N{t') = Y5{t-t') 



(76a) 
(76b) 



with T = 2(j 2 /t [42]. This comparison suggests that 
ujs experiences a "frictional force," of amplitude — (ujs — 
ZUs)/t, that tends to bring that vector back to its origi- 
nal value in spite of the fluctuations. This "force" is pre- 
sumably related to the physical constraint L to t — const., 
although we do not explore the nature of that connection 
here. 

A stochastic realization of an OU process X(t) can be 
generated via [44] 



X (t + At) - X = [X{t) - X] e" At/7 



a 2 (l- e - 2At / T )] V2 n (77) 



where n is a sample value of the unit normal random 
variable, and, in our case, X(t) is one of the compo- 
nents of 0)5. Figure 11 shows an example generated from 
Eq. (77) using the values of {Tj,<7j} extracted from the 
TV-body simulation data of Figure 8. For this example, 
ZJ X and uj y were zero (rotation of the cluster about the 
x-axis) and ZJ Z was set to its initial value. 



We can use these results to rewrite Eqs. (8) and (9) 
in an approximate way that incorporates the effects of 
star-star interactions. Orbits that satisfy the condition 
(52) at t = are assumed to evolve, collisionlessly, in 
response to spin-orbit torques, according to Eq. (9), and 
the contribution of these stars to u>s, which we call u>k, is 
computed as in Eq. (8b). In the case of orbits that do not 
satisfy (52), no attempt is made to follow their detailed 
evolution. Instead, these orbits are assumed to make a 
collective, stochastic contribution to lj s which is modeled 
as an Ornstein-Uhlenbeck time series, [wsjciu, evaluated 
numerically via Eq. (77). The parameters (r, cr) that ap- 
pear in that equation are estimated as described above. 
These two contributions to u>s{t) are then added, and 
the evolution equation for S is written as 



S 



uj s x S, 

U K + Ks]oU- 



(78a) 
(78b) 



Eq. (78) ignores the effects of frame-dragging on stars 
in the "collisional" region. It therefore rules out the 
possibility that differential precession of stars in this re- 
gion could distribute their Lj vectors uniformly about 
S, causing their net torque on the SBH to drop, as oc- 
curs in the collisionless regime (Figure 3). While we can 
not rigorously defend this approximation, we can state 
a more basic set of physical assumptions from which it 
follows. Consider a star whose orbit evolves in response 
both to frame-dragging from the SBH and gravitational 
encounters from other stars. Idealize the encounters as 
occurring at discrete times separated by ~ Trr. Between 
encounters, the line of nodes precesses uniformly at the 
Lense-Thirring rate, by an amount 



AO w uj lt x T RR . 



(79) 



If the effect of an encounter is to randomly select a new 
n - that is, if memory of the previous f2 is completely 
erased after one relaxation time - then the mean change 
in after many encounters will be just Ail. Finally, 
if wltTrr *C 1, then AO <?C 2ir, implying a negligible 
amount of differential precession about S even after ar- 
bitrarily long times. 

A similar argument [45] can be used to derive the drift 
velocity of an electron that is subject to a fixed electric 
field (the SBH torque) and to random collisions (grav- 
itational encounters); the finiteness of the drift velocity 
(nodal angle fl) follows from the assumption that colli- 
sions restore v to a thermal distribution, i.e. that knowl- 
edge of the velocity accumulated prior to the collision 
is lost. As is well known, under some circumstances the 
charged particle retains memory of the velocity it had be- 
fore its collision leading to "persistence-of-velocity" cor- 
rections. We expect our model for spin evolution to be 
similarly limited in its applicability, although we post- 
pone a more thorough understanding of such issues to a 
later paper. 

Some support for this physical picture is provided by 
Figure 2 of Merritt et al. [32], which shows a set of 
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short, numerical integrations of TV-body systems subject 
to frame-dragging torques. The lower left panel in that 
figure corresponds to the case Trr ss w^, and the lower 
right panel to the case Trr. -C i^lt- ^ n the first case, 
stars exhibit a finite amount of nodal precession in spite 
of the encounters, as implied by Eq. (79), while in the 
latter case, encounters appear to remove all traces of a 
net advance of 0. 

Stars near the inner edge of the "collisional" region, 
a > ok, have wlt ^ ^rr: an d for these stars, Afl is not 
necessarily small, as in the lower-left panel of the figure 
just cited. Given that the vaue of aj< is itself uncertain, 
the additional uncertainty due to the evolution of orbits 
in this "transition zone," ax < r <C ol, seems acceptable. 
We note that these uncertainties mimic uncertainties in 
twisted accretion disk models about the location and ra- 
dial extent of the "warp" that determines the torque on 
the SBH, as discussed in more detail in §VIC. 




10* 10' 
M, (r<0.1pc) (M ) 

FIG. 12: SBH precessional periods that would result from 
torquing by stars that orbit outside the rotational influence 
sphere of the SBH, assuming / = 1/2 and x = 1- In the left 
panel, M, = 10 6 Mq is assumed, and the thin and thick lines 
correspond to m* = IMq and m* = 10.M© respectively. The 
right panel assumes the relation (22) between r m and M, and 
m* = IMq. The points were computed from a Monte-Carlo 
model that approximates the observed dependence of 7 on 
M.; dashed lines are for constant 7, as labelled. 

In nuclei with 9k <C 1 (Figure 7), differential preces- 
sion of the stars that contribute to wk will cause their 
torque to die away before the direction of S has changed 
appreciably (Figure 3). Subsequent evolution of S will 
be determined by all the other stars. In a nucleus de- 
scribed by the N(a,e) of Eq. (19), the contribution to 
cjs from those stars (ignoring stochastic fluctuations) is 
given by the integral (39), after restricting the region of 
integration to the complement of (53). The result is 
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3-7 



ok 
r a 



-5/2 



8 y? (3- 7 )(6-7) r( 7 + i) 
17 J 3 2^ (l-2 7 ) 2 r(7-l/2) 



, (80a) 
(80b) 



where {M , r } = {M .i, 0.1 pc} in low-mass galaxies and 
{M ,r } = {2M.,r m } in high-mass galaxies. 

Figure 12 plots this contibution to the SBH preces- 
sional period (that is, to the mean value of 27r/[cjs]ou) 



as a function of nuclear parameters. The results turn out 
to be strongly dependent on 7, the slope of the nuclear 
density profile, so we consider that parameter in more 
detail. Observationally, 7 exhibits a substantial scatter, 
but there is a well-defined mean trend with galaxy lu- 
minosity, at least among the bright galaxies for which 
7 is well-determined [20, 21]: 7 is smaller in the nuclei 
of brighter galaxies. Using standard expressions for the 
mass-to-light ratio of old stellar systems [46] and for the 
mean ratio of SBH mass to galaxy mass [47] , we can write 
this mean relation as 



(7) w 2.0- 1.1 log 



M. 



10 (l<r\W . 



M. > 10'' b M.. (81) 



The right panel of Figure 12 shows a Monte-Carlo distri- 
bution of points generated from this relation, assuming 
a dispersion of 0.25 in 7 at each M,. (Values of 7 < 1/2 
were excluded for the reasons given above.) While the 
scatter is large, there is also a steep trend in the sense of 
smaller precessional periods at lower M,. 



C. Examples of stochastic evolution 

We first consider a dense nucleus in a low-mass galaxy, 
M ga i w 10 9 7W Q : we set M. = lO 6 ^!©, 7 = 2 and 
M0.1 = 2 x 1O 5 A^0. The characteristic radii relating to 
orbital coherence times are given for this nuclear model 
by Eqs. (47) and (55): 



as « 0.27 mpc 



a K 



0.17X 1 



i/2 ( m±_ 



-1/4 



mpc. 



(82a) 
(82b) 



The total angular momentum associated with stars orbit- 
ing close enough to the SBH that frame dragging dom- 
inates self-interactions, a < <Ik(1 — e 2 ) -3 / 4 , is given by 
Eq. (59): 



(83) 



Evidently, an insignificant number stars are in the colli- 
sionless regime. Almost all stars in the collisional regime 
will also have a > as; for these stars, the coherence time 
related to changes in eccentricity (Appendix) is given by 
Eq. (42): 



i coh = i M (e = l/2)«1.0x 10 4 



0.1 pc 



1/2 



yr. (84) 



The incoherent, resonant relaxation time corresponding 
to this i co h, Eq. (A3), is then 



Trr. 



m 



3 x 10' 



1 mpc 



3/2 



yr 



(85) 



which is the time associated with random changes in or- 
bital eccentricities. The 2d resonant relaxation time, in 
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either the coherent or incoherent regimes, is given by 
Eq. (50): 



T 2d RR w 1.1 x 10 4 



-1/2 



( a 



\rnpc 



yr 



(86) 



which is the time scale associated with changes in orbital 
inclinations. As expected, T2dRR < Trr,m- 

The average, spin precessional period of the SBH due 
to torquing by stars in the collisional regime is given by 
Eq. (80): 



2tt 



3.7 x 10 6 



±_ 

0.5 



X 



3/4 



W* 
I Mr. 



-3/8 



yr. (87) 



Even assuming a low degree of net rotation of the cluster 
(/ 0.1), spin precessional periods are predicted to be 
shorter than ~ 10 8 yr. 



x. o 




5x10' 



FIG. 13: x-component of the spin vector % = cS/GM?, de- 
rived as the solution to the stochastic differential equation 
(78), with various values of r: r — 10 3 yr (thinnest, black), 
10 4 yr (red), 10 5 yr (green), 10 6 yr (thickest, blue). Other 
parameters are m* = IMq, f = 0.1, 8o = 50°. 

Since T^^rr ~ t <C 2tt/ujs, we expect that the time 
evolution of S will depend only weakly on the value cho- 
sen for t in Eq. (77), as long as the inequality is main- 
tained. This expectation is confirmed in Figure 13, which 
shows the evolution of the rr-component of S in a set of 
integrations with different r and with / = 0.1, 8q = 50°. 
Only when r is unphysically large, > 10 6 yr, and com- 
parable with the spin precessional period does S(t) show 
an appreciable dependence on it. 



TABLE I: Parameters for Figure 14 
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FIG. 14: x and y components of the spin vector x — cS/GM? 
for a set of integrations of fixed Mo.i = 2 x lO 5 ^© but dif- 
ferent iV and m*: = 0.1.M© (thick/black), m* = 1A4© 
(dash-dotted/blue), and m* = lOA^o (dotted/red). Other 
parameters of the models are given in Table I and in the text. 



Figure 14 shows the evolution of the SBH spin vector 
in a set of integrations with / = 0.1 and with different 
to*; the number of stars was varied in order to keep Mo.i 
fixed at 2 x 1O 5 A^0. The correlation time t was fixed at 
10 5 yr. The number, Nk, and total mass, Mr, of stars 
in the collisionless regime are listed in Table I. As m* 
is increased (i.e. N is decreased), the spin precessional 
period drops, and the dependence of S on time exhibits 
more stochasticity. Both effects are consequences of the 
increasing number of stars in the collisional regime. 

Next we consider the nucleus of an intermediate-mass 
galaxy, M ga i ~ 10 n 7W o . We set M. = 1O 8 A^ and 
777,* = IA4q. The SBH influence radius is r m 35 pc 
(Eq. 22). A typical value for 7 would be - 2 (Eq. 81), 
but given the large scatter in this parameter, we consider 
a range of values. Proceeding as before, we find from 
Eqs. (48) and (56): 



and from Eq. (56) 



«K 



0.21 mpc (7 = 1) 
0.016pc (7 = 2) 



0.16x 2/5 pc (7 = 1) 
0.042x 1/2 pc (7 = 2). 



(88) 
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The angular momentum of stars in the collisionless 
regime is given by Eq. (58) as 



(7 = 1) 



0.20x- 1/4 (^) (7 = 2) 



(90) 



Ok approaches unity for sufficiently large values of / and 
7. As in the previous example, almost all stars in the col- 
lisional regime have a > as, and for these stars Eq. (43) 
gives 

/ \-l/2 

i coh « 3.3 x 10 7 ( yr (7 = 1) 

« 9.3 x 10 5 L^j ' yr (7 = 2). (91) 



Similarly 



Trr„m « 3 x 10 7 



and 
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The spin precessional period due to torquing by stars in 
the collisional regime is 

- « 7.2 xW-/-) V /5 yr (7 = 1) 
w s V - 5 / 



5.0 x 10 



8 ( J_ 

0.5 



X 3/ V (7 = 2). 



(94) 



Note the strong dependence of this time on 7 (Figure 12). 

Figures 15 and 16 show the evolution of S, and of the 
angle between S and L, where L is the angular momen- 
tum of stars in the collisionless regime, for models with 
7 = {1,2}, / = 0.5, 80 = 50°, and various values of 
X- In both sets of model, the long-term precession rate 
of the SBH depends modestly on x, and strongly on 7, 
as expected from the relations (94). There is an initial 
phase in which the stars in the collisionless regime near 
the SBH differentially precess about the nearly-fixed S; 
the length of this phase is ~ 2 Gyr for 7 = 1 and ~ 0.1 
Gyr for 7 = 2. During this time, S reacts somewhat 
to the changing Wk, before settling in to a more regular 
precession (driven by [(*>s]ou) a ^ later times. 

Particularly in the case 7 = 2, it is clear that the angle 
between S and L never reaches zero. In this model, the 
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FIG. 15: The upper panel shows evolution of the x compo- 
nent of the SBH spin in a set of integrations with M, = 
10 s M @ , 7 = 1, / = 1/2 and O = 50°. The different 
curves correspond to different values of the dimensionless spin: 
X = 1.0 (solid/black), x = 0.8 (dash-dotted/red), x = 0.6 
(dashed/green), and x = 0-4 (dotted/blue). The lower panel 
shows the angle between S and L, the total angular momen- 
tum of stars in the "collisionless" regime (i.e. near the SBH). 



time for stars within the rotational influence sphere to 
differentially precess about S is ~ 5 x 10 7 yr, only a few 
times smaller than the SBH precessional period; thus the 
differential precession can never quite "catch up" with 
the changing spin direction. In the case 7 = 1, the ratio 
between these times is more than a factor 10 and the two 
vectors can nearly align. 

Figure 17 plots the evolution of S, L and u>s for 
two models with 7 = 2, x = 1 and two values of 
60 = {70°, 120°}; other parameters are as in Figure 15, 
and L and ujg refer to stars in the collisionless regime 
only. In these models, L is strongly misaligned with S 
initially, and its direction evolves in a very complicated 
way at early times before nearly aligning with S. 



VI. DISCUSSION 

A. Observations of nuclear rotation 

Evolution of SBH spins via the mechanism discussed 
here is a strong function of the degree of rotation of the 
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FIG. 16: Like Figure 15 but for 7 = 2. 
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FIG. 17: Evolution of S, L and ojs in two models with x = 1 
and #o = 70° (upper) and 9q — 120° (lower); other parameters 
are as in Figure 15. The quantities L and lj refer to the 
total angular momentum, and the spin precessional vector, 
due to stars in the collisionless regime. Open/filled circles are 
initial/final values; elapsed time is 0.5 Gyr. 



nucleus, at distances r g <C r <C r m from the SBH, where 
r g is the gravitational radius (Eq. 1) and r m the radius of 
influence (Eq. 21). Observational constraints on the de- 
gree of rotation at such small radii tend to be weak. The 
nucleus of the Milky Way is the closest. Figure 18 plots 
line-of-sight, rotational velocity data for binned samples 



of stars at projected distances < r m 2.5 pc ~ 65" from 
Sgr A* [15, 48]. Also shown for comparison are rota- 
tional velocity curves predicted by the spherical models 
used here (§111); recall that the degree of net rotation 
in those models is set by the parameter /, with / = 1/2 
corresponding to maximal rotation. The Milky Way data 
are consistent with all values of / but the available data 
extend inward only to ~ 0.2r m , well outside the region 
that would contribute most of the torque to a spinning 
SBH. 

The complexity of the velocity data has led to sug- 
gestions [12] that the Milky Way nucleus consists of a 
superposition of different structures with different axes 
of rotation. While the stars contributing to the velocity 
data in Figure 18 are mostly old, two disklike structures 
of young stars - the clockwise disk discussed above, and 
another (the "counter-clockwise disk") [7], both at ~ 0.1 
pc - are known to rotate about axes that are separated 
by ~ 60° and both disks are inclined with respect to the 
large-scale symmetry plane of the Galactic disk. 

Projected Radius (pc) 
0.01 0.1 




0.01 0.1 1 10 

Projected Radius (arc sec) 

FIG. 18: Solid lines are predicted, observed rotational veloc- 
ity curves for a spherical stellar system around a SBH of mass 
4 x W 6 Mq at a distance of 8 kpc, similar to the Milky Way 
nucleus. Curves are labelled by the fraction / of flipped or- 
bits; 7 = 1.5 was assumed. The nucleus is assumed to be 
observed from a point in a plane perpendicular to the axis of 
nuclear rotation, and the velocities have been averaged over a 
perpendicular distance of 1 arc second w 0.04 pc. The dashed 
curve is the Id velocity dispersion in the case that f — 0. Tick 
marks label cik, Eq. (56), and ol, Eq. (27), assuming x = lj 
Mo.i = 10 Mq, and m* = IMq. Circles with error bars are 
measured, binned, line-of-sight mean velocities of stars from 
[48] (open) and [15] (filled). 

The Local Group dwarf galaxy NGC 221 (M32), at 
a distance of ~ 700 kpc [49], also appears to contain a 
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SBH, with mass that is poorly constrained but believed 
to be similar to that of the Milky Way SBH [50]. The 
line-of-sight mean stellar velocity in this galaxy is roughly 
constant with radius, v YOt ~ 60 km s _1 [51, 52], inside a 
projected radius of <~ r m m 3 pc, compared with a line- 
of-sight velocity dispersion a ~ 75 km s _1 , suggesting 
/ < 0.5. The resolution in this case is ~ 0.3 pc <~ 0.1r m . 

Beyond the Local Group, massive galaxies are the best 
prospects for spatially resolving a region smaller than r m 
due to the scaling of r m with galaxy mass. A region of 
size 0.1r m in a galaxy at distance D has angular extent 

«, ( „, lrm) ,0»07( T ^-) O ( T5 £-j)" (95) 

where a s=s 0.56 (Eq. 22). Observations rarely exceed the 
~ 0".l resolution of STIS on the Hubble Space Telescope 
and so little data are available on scales much less than 
0.1r m for galaxies beyond the Local Group. Neverthe- 
less, a number of nearby galaxies exhibit strong nuclear 
rotation, w rot ~ a on the smallest resolvable scales. Some 
examples (NGC number, followed by the radius of the re- 
solved region, expressed as a fraction of r m ) are: NGC 
3115 (0.04) [53], NGC 3377 (0.20) [54], NGC 3379 (0.10) 
[55], NGC 4342 (0.12) [56], NGC 4258 (0.17) [57]. Data 
like these are at least consistent with the presence of sig- 
nificant nuclear rotation on spatial scales <C r m although 
of course they do not compel it. As noted above (§111), in 
giant galaxies, the radius containing an orbital angular 
momentum equal to S is expected to be <~ 0.01r m — 0.1r m 
in the case v rot ~ a. 



B. Model constraints on the degree of nuclear 
rotation 

Given the weak observational constraints on rotation 
of galactic nuclei, it is interesting to ask what various 
models of nuclear evolution predict. 

Rotation arises naturally if stars formed in a thin 
gaseous disk around the SBH [58, 59], before later being 
scattered (say) into more spheroidal structures. Star for- 
mation requires that the gas disk be dense enough for its 
internal gravity to overcome shearing and tidal stresses 
from the SBH. Steady-state accretion disk models [60] 
suggest a minumum radius for star formation of [61] 
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of 10 2 pc is similar to the inner radius of the 



io- m 14/27 f-L) 
V0.03/ Vo.i/ 



8/27 



O.ILp 



-8/27 



/ M. 



\ 10 8 A4n 



1/27 
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where a is the standard viscosity parameter [62], L is 
the luminosity due to gas accretion onto the SBH, e is the 
accretion efficiency defined by L — eMc 2 , and Le a 1.4x 
1O 46 (M./1O 8 7W ) erg s" 1 is the Eddington luminosity. 
The predicted dependence of r min on M, is extremely 
weak. 



young "clockwise disk" of stars at the Galactic center 
[7, 9, 63]. However there is currently no evidence of an 
accretion disk [64, 65] and the low luminosity of Sgr A* 
places strict limits on its current rate of gas accretion 
[66]. Attempts to explain the formation of the young 
stars usually invoke instead the recent infall and tidal 
shearing of a massive gas cloud. Numerical simulations of 
this scenario [67-70] have confirmed that formation of a 
disk from which stars subsequently fragment is possible if 
the initial conditions (cloud mass, density, temperature; 
orbital parameters) are correctly chosen. Star formation 
in these models takes place as close as <~ 0.01 pc to the 
SBH. Given the small number of published simulations, 
and the fact that they were motivated by a desire to 
reproduce the known properties of the stellar disk, is not 
clear whether different initial conditions might allow star 
formation much farther in. 

Figure 2 suggests that for galaxies with M, > 1O 7 A^0, 
a K <; 10~ 2 pc. For these galaxies, restricting the region 
of significant rotation to r > 10~ 2 pc would reduce some- 
what the contribution to ojs from stars in the "collision- 
less" regime but would not change the implied rate of 
steady SBH precession due to stars beyond cxr, as given 
by Eq. (80). In the case of low-mass galaxies, remov- 
ing stars inside ~ 10~ 2 pc would essentially turn off the 
collisionless contribution to S and increase the SBH pre- 
cessional period due to stars in the collisional regime by 
an approximate factor (aK/?*min) 1 / 2 1 ■ 

Another possible source of nuclear rotation is inspiral 
of a massive object, which transfers its orbital angular 
momentum to the stars via dynamical friction before be- 
ing captured by the SBH (say). Assume that the inspi- 
ralling object has a mass to,, where to* <C to, -C M,. 
Assuming a circular orbit, a decrease in orbital radius of 
Ar implies a transfer to the stars of angular momentum 



AL = 



GM. 



Ar. 



(97) 



We want to compare this with the maximum, net angular 
momentum that could be associated with the stars in a 
shell of thickness Ar: 



AL* = f x r x 



GM. 



x 4irr 2 p(r) Ar (98) 



where / < 1 depends on the morphology of the nucleus 
and the distribution of stellar orbits. Thus 
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In terms of the density model adopted here for low-mass 
galaxies, this can be expressed as 



dL 



1 



TO, 



2(3- 7 )/M .i VO.lpc 



7-3 



(100) 



Since 7 < 3, this result implies that the largest fractional 
increase in orbital angular momentum occurs for stars 



19 



0.3 



CD 
> 



0.1 













^ \ 
~$ ■ 






- 








v rot 

. / \ . y 




I 




0.1 1 




10 











FIG. 19: Rotation in the galaxy merger simulations of Gua- 
landris & Merritt (2012) [71]. Solid lines are the one- 
dimensional velocity dispersion, dashed lines are the rotation 
velocity about the z axis (the orbit of merging galaxies was 
the x — y plane). Thick (thin) lines are for a merger on a circu- 
lar (eccentric) orbit. The galaxy mass ratio was 1:3 and each 
galaxy hosted a SBH with M, = 0.005 M ga i axy . The mod- 
els contained four different stellar masses representing an old 
stellar population; stars from all mass groups were weighted 
equally in constructing this figure. The dotted line shows a 
Keplerian rise in velocity dispersion and r m is the radius con- 
taining a stellar mass equal to twice the combined mass of the 
two SBHs. 

nearest the SBH. The model must break down at radii 
where the enclosed stellar mass is less than ~ m., or 

r»0.1p=(j£) ■ HOD 

For example, setting m. = lO 3 ^© (an "intermediate- 
mass black hole"); M .i = lO 5 ^©; and 7 = 2, we find 
fmin ~ 1 mpc. At this radius, \dL/dL±\ is maximized 
and equal to l/[2(3 — 7)/] which can be of order unity. 

Somewhat larger changes in nuclear structure and 
kinematics would result from the dissipationaless (gas- 
free) merger of two galaxies containing comparably- 
massive SBHs [72] . This is a likely model for the origin of 
the cores that are observed at the centers of galaxies with 
M, > 10 7 5 A^o [73]. Orbital motion of the two galaxies 
would imprint rotation on the stars in the merged galaxy, 
but the binary SBH also displaces a mass in stars of order 
its own mass via the gravitational slingshot [74] ; the net 
rotation of the stars left behind depends in a complicated 
way on this process. Figure 19 shows results extracted 
from perhaps the highest-resolution study to date of this 
interaction [71] . The figure shows the velocity dispersion 
and rotation velocity profiles of the merged galaxy at the 
time when two SBHs coalesce. Despite limited statistics 
at small radii, it is clear that such merger products may 
have a noticeable degree of rotation well within the radius 



of influence, corresponding roughly to / ~ 0.1 — 0.2. 

Both infall of gas clouds and inspiral of massive com- 
pact objects could occur episodically. In the case of spiral 
galaxies, the example of the Milky Way with its young 
stellar disks suggests that accretion events might occur 
roughly once per 0.1 Gyr [75]. In the case of massive 
galaxies, which tend to be gas-poor, large-scale simula- 
tions of dark-matter clustering suggest that the mean 
time between galaxy mergers varies from ~ 0.2 Gyr at 
a redshift z = 10 to ~ 10 Gyr at z = 1 with a weak 
dependence on galaxy (i.e. dark halo) mass [76]. As- 
suming that all or most galaxies contain nuclear SBHs, 
this would also be roughly the time between insertion of 
secondary SBHs into the nucleus [3]. These times are 
comparable with the time scales for spin precession de- 
rived here (e.g. Figure 12), suggesting that the evolution 
of SBH spins due to frame dragging may also be episodic 
in nature. 

Both sorts of infall event are likely to occur from essen- 
tially random directions, so that the increase over time of 
the net rotation of the nucleus will have the form of a ran- 
dom walk. Futhermore, both sorts of event can change 
the magnitude of S: if accretion of the gas by the SBH 
occurs [61]; or if the inspiralling body coalesces with the 
SBH [77]. 



C. Comparison with accretion disk torquing 

Interaction of a spinning SBH with a misaligned, 
gaseous accretion disk is driven by the same frame- 
dragging torques modelled here. The accretion disk prob- 
lem has been extensively studied [2, 5, 78-81], in part 
because the radio jets that power the classic, double ra- 
dio sources are believed to be launched perpendicularly 
to the inner accretion disk. Both the long-term (<~ 10 8 
yr) stability of jet directions in some active galaxies, as 
well as the evidence for jet precession in others, is prob- 
ably linked in fundamental ways to SBH - accretion disk 
interactions. 

Here we sketch the points of similarity and difference 
between spin evolution driven by a misaligned accretion 
disk and by a rotating stellar nucleus. We emphasize 
that the latter case is generic - SBHs appear always to 
be embedded in stellar nuclei - while nuclear activity, 
hence accretion disks, exist in only a small subpopulation 
of galaxies. The fact that accretion disk torques have 
received essentially all the attention until now is probably 
a consequence of the easy observability of the jets. 

Given an assumed structure for the disk (surface den- 
sity and inclination as functions of radius), the instanta- 
neous evolution equation for S is essentially Eq. (8), after 
setting orbital eccentricities in that equation to zero and 
identifying Lj with the angular momentum of a discrete 
element of gas. Such models typically assume a disk that 
is thin and initially planar. Differential precession then 
ensues near the SBH; at radii r -C - defined, as in 
Eq. (27), as the radius containing an angular momentum 
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equal to S - the gas precession time is short compared 
with that of the SBH. 

The value of tl can be computed given a model for the 
disk surface density. For instance, the steady-state disk 
models referred to above [60] imply [81] 
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where a, e and Le are defined as in Eq. (96). In active 
galaxies, all quantities in parentheses aside from the fac- 
tor containing M, are of order unity and so 10~ 25 pc < 
tl < 10 -0 5 pc. This is somewhat smaller than the value 
of ol as plotted in Figure 2, at least in massive galax- 
ies; in other words: accretion disks, when present, are 
likely to dominate the angular momentum distribution 
near the SBH, justifying the neglect of stellar torques in 
these galaxies. In smaller galaxies, Figure 2 suggests that 

Differential precession causes gas near the SBH to at- 
tain a mean L that is aligned with S, as in the stel- 
lar case, but gaseous viscosity also ensures that the gas 
returns to a thin disk, coincident with the SBH equa- 
torial plane (the "Bardeen-Petterson effect" [2]). This 
thin, aligned disk extends outward, not to (where the 
Lense-Thirring time is likely to be very long anyway) , but 
rather to the smaller radius r warp , the radius at which the 
disk plane transitions to its large-radius orientation. The 
warp radius is determined by the condition that the time 
scale for angular momentum diffusion through the disk 
is equal to the Lense-Thirring precession time: 
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Here v is the kinematic viscosity, which also determines 



the accretion rate. The value of ■ 



warp 



is strongly model- 



dependent and still rather uncertain; early estimates (e.g. 
[2, 79]) set r warp w n,, but more recent estimates (e.g. 
[5, 82]) find r warp « (10 2 - 10 3 ) x r g < r L . 

Once alignment of the gas inside r war p has occurred, 
precession of the SBH is driven by gas at r > r warp . An 
expression that is often given for the steady-state SBH 
precession frequency (e.g. [5, 79]) is 
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where L(r < r warp ) is the angular momentum of disk 
gas inside r warp . Uncertainties about the value of r warp 
translate via this expression into uncertainties about the 
precession rate. Equation (104) is similar to Eq. (14) 
for the mutual precession of a SBH and a ring of mat- 
ter, especially when it is recognized that J w S in the 
accretion-disk case. This is at first sight surprising, since 
Eq. (104) appears to relate the precession of the SBH 
to the angular momentum of gas all of which, by as- 
sumption, is fully aligned with the SBH! The justification 



(e.g. [79]) consists of noting that L x wlt oc L(r)/r 3 is a 
steeply falling function of radius, hence only matter near 
the warp is relevant. But this argument underscores the 
very approximate nature of Eq. (104). 

The warp radius plays approximately the same role as 
the radius ok in the stellar case, Eq. (53). The SBH 
precession frequency, Eq. (104) in the gaseous case, be- 
comes Eq. (80) in the stellar case. In the stellar case, 
ok <C gil (Figure 2), just as r warp <C tl (at least if the 
most recent estimates of r warp are correct). 

The continued deposition of matter from a fixed outer 
plane must ultimately align S with the outer L. In many 
models [5, 80, 83], the time scale for this alignment is 
similar to the warp-driven precession time of the SBH, 
i.e. the inverse of Eq. (104). Typical values quoted for 
^aiign lie in the range 10 7 — 10 8 years and it has been ar- 
gued that this alignment is responsible for the long-term 
(10 8 — 10 9 yr) stability of jet directions in many active 
galaxies. Interestingly, we found that complete alignment 
was possible also in the stellar-dynamical case (Figure 4) ; 
differential precession is sufficient to achieve this, even in 
the absence of viscosity. However we argued that a more 
generic outcome in the stellar case is steady precession 
of the SBH, particularly when stellar interactions are al- 
lowed. 

The evolution of S due to the combined influence of a 
misaligned accretion disk and stars is beyond the scope of 
this paper, but we include a few speculative remarks [97] 
Feeding of active galaxies is probably episodic [84, 85]. 
When much, but not all, of the infalling gas has been 
consumed, there may come a time when the precession 
rates due to gas and stellar torquing are comparable. If 
the SBH is still active at this time, accretion-disk-related 
jets should begin to precess roughly in the manner dis- 
cussed here, even if the SBH had previously reached a 
steady-state alignment with the gas. Prolonged, steady 
precession of radio sources might be explained in this way 
[86, 87]. After the gas has been fully consumed, the SBH 
spin can continue to evolve in response to the stars. If the 
gas has been accreted all the way to the event horizon, 
both the magnitude and direction of S will have been 
changed by the gas. 



D. Slowly-rotating nuclei 

Even in a nucleus with negligible net rotation, there 
will still be a nonzero torque on the SBH due to im- 
perfect cancellation of the Lj from the finite number of 
stars. This is obvious, for instance, from Figure 8; the 
components of u>s perpendicular to the mean rotation 
axis of the cluster are zero on average but fluctuate as 
orbits change their Lj due to encounters. Eq. (73) is an 
estimate of the size of those fluctuations and can equally 
well be interpreted as the expected value of uis in a non- 
rotating, isotropic cluster with known N(a). 

It is interesting to ask how large the steady rotation of 
a nucleus needs to be if the net torque exerted on the SBH 
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is to exceed this (fluctuating) value. We estimate the 
torque in a nonrotating nucleus by setting r p m i n = ax 
in Eq. (73); in other words, we conservatively ignore the 
torque from stars within the sphere of rotational influence 
given that they may have differentially precessed about 
S. The result is 
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where {M , r } — {Mo.i, 0.1 pc} in low-mass galaxies and 
{Mo, ro} = {2M # , r m } in high-mass galaxies. Comparing 
a to ujq as given by Eq. (80), we find for the critical 
degree of rotation 
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For instance, in a low-mass galaxy with 7 = 2, 
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For this value of /, the instantaneous time over which S 
changes is 
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We emphasize that this is not a precessional period since, 
by assumption, finite-iV effects dominate and the axis 
about which S is precessing will itself change, in a time 
of order T 2d RR. 

When observing changes in S over very short time 
scales, e.g. human lifetimes, there will also be a time- 
dependent contribution due to the motion of stars along 
their (unperturbed) orbits. This contribution has been 
ignored up till now due to the orbit-averaging of Eqs. 
(8)-(9). 

Figure 20 illustrates the complexity of the evolution 
of S in the case that all the torque on the SBH is due 
to these finite- TV effects. The figure is based on a direct 
iV-body integration of a cluster of 50 "stars," of mass 
5OA^0 each, around a SBH of mass 10 6 A^o an d X = 1- 
Additional details about the initial models are given in 
Merritt et al. (2011) [37]. Over the ~ 2 Myr time span 
of the integration, the SBH spin axis wobbles by about 
one degree. 



E. Experimental determination of black hole spins 

Several authors [58, 88, 89] have suggested that it 
may be possible to infer the magnitude and direction 
of SBH spins from the precession of the angular momen- 
tum vectors of individual stars. Only stars with semima- 
jor axes a < ok are suitable for this purpose, otherwise 
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FIG. 20: Evolution of the dimensionless spin in a direct N- 
body integration [37] . The top panel shows the x and y com- 
ponents of x which is initially parallel to the z-axis. Open 
and filled circles show initial and final orientations. In the 
bottom panel, the dashed (blue) curvs is % x 5 

the dotted (red) 
curve is Xvi an d the solid (black) curve is Xz- 



the changes of Lj due to collisional effects will supersede 
changes due to frame dragging [32, 90]. Since for most 
galaxies ok <C ol (Figure 2), this also means that the 
Lense-Thirring precession times for these stars are much 
shorter than SBH spin precession times, or Lij 3> Uis- 

Assuming that S precesses steadily about the (fixed) 
axis ws, it is convenient to consider the evolution of Lj 
in the reference frame which rotates with the precession 
frequency u>s about the axis of SBH precession, so that 
S is stationary. In this rotating frame, the equation of 
motion for the orbital angular momenta (9a) reads 



OJi = U3j - LJ S 



(109) 



where X/ is the angular momentum in the rotating frame. 
In other words: the Lense-Thirring precession is occuring 
about an axis that does not coincide with the instanta- 
neous direction of S. However, the first derivative of Lj 
in the inertial frame coincides with the value obtained 
without taking into account SBH precession, since it is 
determined by the instantaneous value and direction of 
S. It is only for the second derivative of Lj that the 
difference starts to matter. Given that |u>'- — u)j \ <C \uJj\ 
because we are considering stars that themselves precess 
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much faster than the SBH, it seems unlikely that these 
effects may be detectable in the near future. 

In a nucleus that is sufficiently old, differential pre- 
cession of stars with a < ax will have caused their angu- 
lar momentum vectors to distribute themselves uniformly 
about S, as shown above in several numerical examples. 
This suggests a way of measuring the instantaneous di- 
rection of S via the mean direction of the Lj within ax- 
In the Galactic center, ax ~ 10~ 3 pc, and there are as 
yet no stars with determined orbits in this region. How- 
ever orbital periods for stars with a = 1 mpc are about 
one year, and it is possible that determination of the or- 
bital elements of a few such stars might be feasible over 
a shorter time interval than is required for measuring 
changes in the Lj due to frame dragging. 

VII. CONCLUSIONS 

1. In a galactic nucleus containing a spinning supermas- 
sive black hole (SBH), frame dragging results in mutual 
torques between the stellar orbits and the SBH. The re- 
sult is precession of both the SBH spin, S, and the an- 
gular momentum vectors, Lj, of the individual stellar 
orbits, with S + L to t — S + J2j Lj conserved. For stars 
at a single distance from the SBH, the controlling param- 
eter is the ratio between S and L to t ■ If S 3> L to t , stellar 
orbits precess about the nearly fixed S with the Lcnsc- 
Thirring period; while if L to t ^> S, S precesses about the 
nearly fixed L to t with a period that is shorter by a factor 
S/L tot . The inner parsec of the Milky Way is known to 
contain stellar subsystems having L tot ~ S. 

2. Ignoring interactions between the stars, solutions of 
the coupled equations for S and ij=i,...,jv in spherical 
nuclei reveal two evolutionary modes in the case that 
Ltot > S: continued precession of S about L tot ; or 
damped precession, in which S and L tot come into nearly 
complete alignment after one precessional period of the 
SBH. Even in the first mode, differential precession of 
orbits near the SBH causes their net angular momentum 
to align with S, reducing the torque that they exert on 
the SBH. Subsequent precession of the SBH is driven by 
torques from stars at r > tl, where Tl is the radius en- 
closing a net angular momentum equal to S. 

3. Newtonian interactions between stars can change their 



orbital angular momenta in a time shorter than Lense- 
Thirring processional times. We define the "radius of ro- 
tational influence," ax, around a Kerr SBH as the radius 
inside of which torques due to frame dragging act more 
quickly than torques from the other stars. Typical values 
for this radius are ~ 10 -3 parsecs in dense nuclei like that 
of the Milky Way, increasing to <~ 10° — I0 1 parsecs in 
nuclei containing the most massive SBHs. The angular 
momentum associated with stars in this "collisionless" 
region near the SBH is likely to be much smaller than S 
in nuclei of the smallest galaxies but may be comparable 
to S in massive galaxies. 

4. Interaction between stars at r > ax leaves the total 
angular momentum of these stars unchanged, but results 
in random fluctuations of the individual Lj and hence in 
the torque which they exert on the SBH. We develop a 
stochastic model, based on the Ornstein-Uhlenbeck equa- 
tion, for the torque exerted by these stars and verify it 
by comparison with high-accuracy iV-body simulations. 
We argue that dS /dt can be approximated as the sum of 
two terms: deterministic torques exerted by stars inside 
ax, whose angular momenta evolve solely in response to 
frame-dragging; and a stochastically-fluctuating torque 
due to stars outside «k- 

5. Examples of stochastic evolution of S are presented 
for various nuclear models. Typical evolution consists of 
sustained precession, with periods that are highly depen- 
dent on nuclear parameters, but which are expected to in- 
crease with increasing M.: likely periods are <~ 10 7 — 10 8 
yr for low-mass SBHs in dense nuclei, ~ 10 8 — I0 10 yr for 
SBH with masses - 1O 8 A4 , and - 10 10 - 10 11 yr for 
the most massive SBHs. 
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I 

Appendix A: Time scale for change in orbital eccentricity 

Here we present approximate expressions [16] for the time scales associated with changes in orbital eccentricity due 
to resonant relaxation [33] and evaluate them for the power-law density model used in the text. 

Define the "apsidal coherence time" t co h as the shorter of the two precession times t$ and *m defined in §V, each 
evaluated at typical values of e; say, e s=s 1/2. Comparison of Eqs. (41) and (50) shows that tm s=s T^rr/VN where N 
is the number of stars at r < a: apsidal precession due to the spherically-distributed mass acts more rapidly than ^/N 
torques at changing orbital orientations. For elapsed times short compared with i co h, the torque due to all the local 
stars is therefore nearly constant, and the angular momentum of a typical star responds by changing approximately 
linearly with time. In this "coherent resonant relaxation" regime, all the components of Lj, i.e. both the orientation 
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angles and the eccentricity ej, change with characteristic time TRR ;CO h given by 

P M. 1 ,( a \ 3/2 / M. Y 1/2 f M. \fNY 1/2 /4 x 

This is the same expression as Eq. (50) for T 2 dRR, reflecting the fact that in the coherent regime, both the direction 
and the magnitude of L change on roughly the same time scale. 

For time intervals longer than t co h, the direction of the net field-star torque changes, and evolution of the Lj in 
response to the torques is better described as a random walk. The time scale associated with this "incoherent resonant 
relaxation" is 
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where L c = y/GM t a is the angular momentum of a circular orbit of semimajor axis a, and AL co h is the change in L 
during At = i co h- Setting t co h = tu (he. a > as), this becomes 
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with Cs again of order unity. Eqs. (A3) and (A4) are the appropriate time scales to associate with changes in orbital 
eccentricity in the incoherent regime. 

In the case of 2d resonant relaxation, the relevant coherence time is that associated with changes of the orbital 
planes, i.e. T2dRR- Since AL co h w L c in this case, Eq. (A2) implies that the coherent and incoherent relaxation times 
are approximately the same: no new time scale arises in the incoherent regime for 2d resonant relaxation. 

Comparing the incoherent relaxation times associated with changes in the orientation and magnitude of L, respec- 
tively, we find 

Tmrr 1 , . . x 
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The first of these ratios is manifestly smaller than unity at all radii. The second is only relevant at a < as, i.e. for 

^V<^4= (A6) 



which implies 



again less than unity. On the basis of these inequalities, it is reasonable to equate the correlation time associated 
with fluctuations in ujg with the shortest of the time scales, T*2dRR, as was done in jjVB. 
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